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Technical  Report  Summary 


The  research  of  O.W.  Nuttli  on  regional  wares  can  be  summarized  a::: 

Attenuation.  Attenuation  coefficient  values  for  1-sec  period  Lg 
waves  have  been  determined  for  paths  from  NTS  to  selected  stations 
in  the  western  United  States,  for  paths  from  the  Sahara  test  site 
to  selected  HVSSn  stations  in  Africa,  for  paths  from  the  Bast 
Kazakh  test  sites  to  selected  WWSSN  stations  in  Scandinavia  and 
southern  Asia,  and  for  paths  across  Iran.  Attenuation  also  was 
determined  for  1-sec  period  Pn,  Pg  and  Sn  waves,  and  for  3-sec 
period  Lg  waves,  across  Iran.  Each  of  these  waves  have  been  used 
to  derive  regional  magnitude  formulas  for  m^. 

nt serial nation.  The  discrimination  ability  of  l-sec  period  Lg-wave 
amplitudes  is  low.  For  the  same  m^P)  value,  the  earthquake  Lg 
waves  are  on  average  only  25 %  larger  (0.1  magnitude  unit)  than  the 
explosion  Lg  waves.  However,  Lg-wave  amplitudes  from  earthquakes 
vary  appreciably  with  azimuth,  whereas  Lg-wave  amplitudes  from 
explosions  show  a  near-circular  radiation  pattern.  Therefore,  the 
standard  deviation  of  mfa(Lg)  should  be  larger  for  earthquakes  than 
explosions,  which  may  serve  as  a  regional  discriminant. 

La  Excitation  and  P-Wave  Magnitude  Bias.  If  the  source  medium  is 
taken  into  acoount,  the  excitation  of  1-sec  period  Lg  waves  is  the 
same  for  explosions  of  similar  yield  in  various  geographic  regions. 
That  is,  there  is  no  magnitude  bias  for  m^CLg)  as  there  is  for 
m^(P).  The  Lg  data  can  be  used  to  estimate  nt^(P)  bias,  if 
m^Cteleselsmic  P)  values  are  available.  Taking  eastern  North  Amer- 
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ica  as  a  datum  (assume  m^Lg)  a  m^P)  for  that  region),  the  m^P) 
bias  Is  0.0  for  East  Kazakh  and  0.4  for  NTS,  the  Sahara  and  Iran. 

Regional  Estimates  of  Explosion  Yield.  The  announced  yields  of 
explosions  at  NTS  are  used  to  empirically  determine  a  relation 
between  m^(Lg)  and  log^  yield.  For  hard-rock  sources,  the  empiri¬ 
cal  equation  is  a  quadratic  log  yield,  with  a  slope  of  near  unity 
at  small  yields  (10  kt)  and  of  about  0.6  at  large  yields  (1000  kt). 
For  alluvium  sources  at  NTS,  the  curve  is  almost  linear,  with  a 
slope  of  approximately  0.7  for  the  range  of  1  to  100  kt. 

The  NTS  calibration  curve  for  Lg  has  been  used  to  estimate  yields 
of  selected  explosions  at  East  Kazakh.  From  April,  1976  through 
September,  1980,  the  largest  yield  is  found  to  be  that  of  14  Sep¬ 
tember  1980,  with  an  estimated  yield  of  195  kt,  based  on  Lg  ampli¬ 
tudes. 


The  research  of  B.J.  Mitchell  was  primarily  concerned  with  the 
attenuation  of  waves  at  intermediate  periods  (2-50  s),  and  the  relation¬ 
ship  of  the  attenuation  of  those  waves  with  that  of  regional  phases  at 
higher  frequencies.  Methods  were  determined  by  which  fundamental-mode 
and  higher-mode  Rayleigh  wave  speotra  could  be  utilized  to  determine  Q_ 

P 

as  a  function  of  the  depth.  The  methods  are  important  in  regional  stu¬ 
dies  of  attenuation  because  they  can  be  used  with  a  single  station  and 
single  event.  They  were  used  to  determine  Q  structure  in  the  eastern 
United  States,  the  Colorado  Plateau,  and  the  Basin  and  Range  province  in 
North  America,  and  in  a  region  of  Iran  and  Turkey  and  in  the  Barents 
shelf.  The  results  indicate  that  Q  varies  from  region  to  region  in  the 
upper  orust. 
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If  the  models  obtained  using  Intermediate-period  waves  are  used  to 
predict  the  attenuation  of  higher  frequency  regional  phases  such  as  Lg, 
it  is  found  that  the  models  predict  muoh  lower  Q  (greater  attenuation) 
than  is  observed  for  those  phases.  In  order  to  explain  that  result  it 
is  required  that  the  dissipative  properties  of  the  crust  be  frequency- 
dependent,  higher  frequencies  being  characterized  by  higher  Q  values. 
This  result  is  confirmed  by  both  frequency-domain  studies  and  synthetic 
seismogram  analyses  of  1  Hz  Lg  waves. 
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Discrimination  Potential  of  1-sec  Period  Lg  Waves 
Otto  W.  Nuttli 

The  ratio  of  m^  (obtained  from  teleseismic  P-wave  amplitudes)  to 
(obtained  from  20-sec  period  surface-wave  amplitudes)  has  been  proven 
over  the  years  to  be  a  useful  discriminant  between  underground  explo¬ 
sions  and  earthquakes.  Three  principal  problem  areas  concerning  the 
effectiveness  or  usefulness  of  the  discriminant  are:  1)  For  small  explo¬ 
sions  the  amplitude  of  20  sec-period  surface  waves  is  too  small  to  be 
detected.  This  occurs  at  an  Mg  value  of  about  3,  corresponding  to  a 
contained  explosion  of  approximately  10  kt  in  hard  rock.  2)  Some  large 
explosions  release  tectonic  strain  and  thus  have  an  earthquake-type 
radiation  which  is  superposed  on  the  radiation  of  waves  from  the  explo¬ 
sion.  3)  Some  earthquakes,  particularly  those  in  mid-plate  environ¬ 
ments,  are  associated  with  high  stress  drop.  Thus  their  source  spectrum 
is  similar  in  shape  to  that  of  explosions,  and  their  :  Ms  values  are 
more  explosion-like  than  those  of  earthquakes.  A  few  Asian  earthquakes 
of  this  type  were  described  by  Nuttli  and  Kim  (1975).  The  high  stress 
drop  of  the  New  Brunswick  earthquake  of  9  January  1 982  had  :  M<,  values 

which  fall  in  the  explosion  population  in  Figs.  10.7b  and  10.7c  of  Dahl- 
man  and  Israelson  (1977).  The  mainshock,  which  occurred  at  00hl4m,  had 
an  =  5.7  and  Mg  =  4.8,  and  an  aftershock  on  the  same  day  at  12h53m 
had  m^  =  5.1,  Mg  =  3.8  (O.S.  Geological  Survey,  1982). 

The  second  problem  ordinarily  is  associated  with  explosions  of 
large  yield,  so  that  other  means  can  be  used  to  distinguish  them  from 
earthquakes.  However,  the  tectonic  release  tends  to  distort  yield  esti¬ 
mates  based  on  surface-wave  magnitude. 
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The  third  problem  is  not  so  tractable,  and  may  be  common  for  mid¬ 
plate  earthquakes.  This  matter  is  considered  further  in  the  section  of 
this  report  entitled  "Magnitude  Relations  and  Spectral  Scaling.” 

Because  Lg  is  a  superposition  of  higher-mode  surface  waves,  a 
number  of  investigators  have  considered  it  to  be  potentially  useful, 
when  compared  with  P-wave  amplitudes,  to  discriminate  between  explosions 
and  earthquakes  for  small  events  for  which  20-sec  period  surface  waves 
are  too  small  to  be  detected.  The  reasoning  in  support  of  this  deduc¬ 
tion  is  that  explosions,  which  are  essentially  compressive  radiation 
sources,  will  not  generate  short- period  surface  waves  as  well  as  earth¬ 
quakes  of  the  same  P-wave  magnitude.  Blandford  (1981)  concluded  that 

"the  amplitude  ratio  of  the  maximum  motion  before  S  (P  )  to  the  max- 

n  max 

imurn  after  SQ  (Lg)  is  a  good  discriminant  between  earthquakes  and  explo¬ 
sions  at  regional  distances.  Experimental  data  shows  that  the  ratio  is 
not  seriously  affected  by  site  or  source  geology  or  by  event  depth. 
However,  propagation  effects  can  be  severe,  especially  from  the  USSR  to 
the  South,  and  can  lead  to  the  requirement  for  regionalization."  In 
Figure  33  of  Blandford  (1981)  vertical-component  Lg  amplitudes,  measured 
at  NORSAR  and  equalized  according  to  the  P-wave  magnitudes  of  the  event, 
are  plotted  versus  epicentral  distance  along  with  attenuation  curves 
obtained  for  the  eastern  United  States.  The  Lg  waves  whose  amplitudes 
are  plotted  were  produced  by  nine  Soviet  explosions,  eight  of  which  were 
PNE's  in  the  western  part  of  the  country.  The  remaining  event  was 
located  at  the  East  Kazakh  test  site.  These  nine  data  points  define  an 
explosion  attenuation  curve.  Lg  amplitudes  at  NORSAR  of  two  Austrian 
earthquakes  also  are  plotted  in  the  figure.  The  one,  for  an  earthquake 
of  17  June  1972,  falls  below  the  explosion  curve,  which  puts  it  in  the 
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"explosion"  population.  The  other  data  point,  for  the  earthquake  of  16 
April  1972,  lies  about  0.5  magnitude  unit  above  the  curve,  in  the 
"earthquake"  population.  Thus,  of  the  two  earthquakes  considered,  one 
passes  and  one  fails  the  discriminant  test.  Figure  42  of  Blandford 
(1981)  shows  that  the  Semipalatinsk  earthquake  of  20  March  1976,  of  = 
5.1,  had  an  Lg^max/Pg^max  ratio  about  ten  times  larger  at  MSH  (Iran) 
than  a  15  January  1976  explosion  of  =  5.2. 

Pomeroy  (Rondout,  1983)  reviewed  the  research  on  discrimination 
using  regional  events.  He  noted  that  the  Teledyne  Geotech  group,  con¬ 
sisting  of  Blandford  and  his  co-workers,  found  the  PmaT/Lg  amp  ,’de 

ratio  to  be  a  good  discriminant,  both  in  Asia  and  the  United  STatt  On 

•» 

the  other  hand,  Pomeroy  noted  that  Bennett  and  Murphy  of  SJ  found  Vi¬ 
olent  scatter  in  the  amplitude  ratio  data  for  western  United  ^.^ces 
earthquakes  that  the  error  bounds  of  the  earthquake  data  enveloped  the 
explosion  data  of  NTS  events,  and  concluded  that  the  Pn/Lg  amplitude 
ratio  was  not  a  promising  discriminant.  In  review  of  Mykkeltveit  and 
Husebye's  work  at  NORSAR,  Pomeroy  noted  that  earthquakes  and  explosions 
in  the  western  and  central  portions  of  the  USSR  showed  large  scatter  in 
the  P/Lg-type  amplitude  ratios,  and  the  authors  could  not  find  any 
disoriminant  of  this  type  that  could  be  applied  to  their  data  with  any 
degree  of  success.  Commenting  on  his  own  research,  Pomeroy  observed 
that  western  United  States  earthquakes  and  NTS  explosions  as  recorded  at 
the  Tellowknlfe  array  in  Canada,  do  separate,  but  with  considerable 
overlap,  when  the  ratio  of  the  Lg  amplitude  to  the  first  significant 
amplitude  in  the  P  wave  train  is  utilized.  Pomeroy  and  Nowak  (1978a,  b) 
noted  that  the  Lg/P  amplitude  ratio  data  for  the  explosion  SALMON  in  the 
eastern  United  States  fell  below  similar  earthquake  data,  but  that  for 
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similar  explosions  in  the  western  Soviet  Union,  the  ratios  for  earth¬ 
quakes  and  explosions  were  the  same. 

Research  done  under  Contract  F49620-79-C-0025  that  is  relevant  to 
the  regional  discrimination  problem  was  reported  by  Nuttli  (1981).  This 
study  considered  16  explosions  and  6  earthquakes  in  the  western  and  cen¬ 
tral  Asia  portion  of  the  USSR.  Short-period  P-wave  amplitudes  were 
measured  (as  done  conventionally  for  determining  m^),  a^on8  with  1-sec 
period,  vertical-component  Lg  amplitudes  from  WWSSN  seismograms  at  13 
stations  in  Germany,  Scandinavia  and  southern  Asia. 

The  variation  of  P-wave  amplitude  with  epicentral  distance,  after 
all  events  were  equalized  to  m,,  =  5.0,  showed  a  scatter  of  approximately 
ten  times  (one  m^  unit)  over  the  distance  range  of  about  5°  to  30°. 
Scatter  of  the  explosion  amplitudes  was  comparable  to  that  of  the  earth¬ 
quake  amplitudes.  The  Veith-Clawson  (1972)  amplitude  attenuation  curve 
provided  a  good  average  fit  to  the  data  over  the  distance  range  con¬ 
sidered. 

The  Lg  amplitudes  individually  had  to  be  corrected  for  the  propaga¬ 
tion  path  from  the  source  to  the  station,  because  the  coefficient  of 
anelastic  attenuation  is  different  for  each  path,  and  because  small 
errors  in  the  assumed  value  of  this  coefficient  can  lead  to  orders-of- 
magnltude  errors  in  calculated  Lg  amplitudes  at  large  regional  dis¬ 
tances.  In  this  study  the  attenuation  coefficient  of  1-sec  period  Lg 
waves  was  determined  from  the  variation  of  the  wave  frequency  in  the  Lg 
coda  with  travel  time  (Herrmann,  1980),  as  described  in  the  section  of 
this  report  dealing  with  Lg  attenuation  in  Eurasia. 
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By  use  of  the  Veith-Clawson  (1972)  P-wave  amplitude  attenuation 
curve  and  of  observed  Lg  attenuation  coefficient  values  for  each  of  the 
90  source-to-station  paths,  theoretical  values  of  log  (A^/Ap)  were  oal- 
culated  for  the  90  path  lengths,  assuming  that  all  events  were  earth¬ 
quakes.  These  theoretical  values  were  compared  with  the  observed  values 

of  log  (A,  /A  ) .  The  deviation,  or  difference  between  the  log  (A.  /A  ) 
°  — g  P  Lg  p 

of  observed  values  and  the  log  (A^g/Ap)  of  calculated  values,  is  plotted 
in  Fig.  1.  From  the  figure  it  can  be  seen  that,  in  general,  the  devia¬ 
tion  for  the  explosions  is  negative  and  for  the  earthquakes  is  near  zero 
or  positive.  This  indicates  that,  on  the  whole,  the  ratio  of  A^  /A  is 
larger  for  earthquakes  than  for  explosions.  However,  the  scatter  and 
intermixing  of  explosion  and  earthquake  data  points  in  the  figure  is 

enough  to  conclude  that  the  amplitude  ratio  A.  /A  of  1-sec  period  waves 

Lg  p 

recorded  at  regional  distances  by  WHSSN  vertical-component  seismographs 
is  not  an  effective  discriminant.  That  is,  too  many  of  the  events  could 
be  mlsldentifled. 

In  the  process  of  compiling  information  for  this  report  we  have 
observed,  on  the  basis  of  a  few  events  only,  that  Lg  amplitudes  from 
earthquakes  show  greater  variation  with  azimuth  than  do  those  from 
explosions  located  in  the  same  region  and  recorded  at  the  same  stations. 
This  suggests  that  Lg  has  a  circular  radiation  pattern  for  explosions 
and  a  lobed  pattern  for  earthquakes.  If  this  observation  proves  true, 
the  standard  deviation  of  a^CLg)  of  explosions  will  be  less  than  of 
earthquakes.  This  idea  will  be  pursued  in  further  studies,  to  determine 
if  it  can  serve  as  an  effective  discriminant. 

On  the  basis  of  the  present  research,  however,  we  conclude  that  the 


lOQ  (A  Lg  A  p)observ«d  “ 
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Fig.  1.  Difference  between  Obeerved  and  Predicted  Values  of  log  (A^g/Ap)  for 
Earthquakes  and  Explosions  in  Eurasia.  If  Lg  waves  frost  earthquakes 
are  not  excited  as  strongly  as  for  earthquakes,  the  explosion  data 


points  should  have  ordinates  of  less  than  zero. 
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scatter  of  A^/Ap  values  for  explosions  and  earthquakes  in  the  western 
and  central  Asia  portions  of  the  USSR  Is  too  large  to  make  the  ratio 
useful  as  a  discriminant.  This  conclusion  is  in  agreement  with  that  of 
most  other  investigations,  except  for  Blandford  and  his  colleagues  at 
Teledyne  Geotech. 
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Magnitude  Relations  and  Spectral  Scaling 
Otto  W.  Nuttli 


Magnitudes  play  an  Important  role  in  discrimination  and  yield  esti¬ 
mation,  both  in  the  application  of  the  m^Mg  discriminant  and  in  the  use 
of  m^P),  m^(Lg)  and/or  Mg  to  estimate  yields.  Therefore,  it  is  of 
value  to  understand  the  factors  which  affect  the  estimates  of  the  magni¬ 
tudes  and  of  the  inter-relations  among  them. 

f 

Research  performed  under  Contract  F49620-79-C-0025  considered  vari¬ 
ous  aspects  of  these  problems,  including  the  relation  of  ^(P)  to  M„  for 
mid-plate  and  plate-margin  earthquakes  (Nuttli,  1983a,  b),  the  relation 
of  m^Lg)  to  Ml  (Herrmann  and  Nuttli,  1982),  and  the  so-called  magnitude 
bias  of  m^P)  caused  by  anomalous  attenuation  of  teleseismic  P-wave 
amplitudes  in  the  asthenosphere  beneath  the  source  (Nuttli,  1982a,  b). 

For  more  than  ten  years  it  has  been  known  that  the  versus  Mg 
curve  used  to  separate  the  explosion  and  earthquake  populations  varies 
regionally.  In  general,  plate-margin  earthquakes  have  larger  Mg  values 
than  mid-plate  earthquakes,  if  both  have  the  same  mb  value  (even  if  P- 
wave  magnitude  bias  is  taken  into  account).  The  n^^g  relations  were 
determined  empirically,  and  varied  somewhat  with  author  and  with  type  of 
instrument  used  to  obtain  the  magnitudes. 

Nuttli  (1983a)  used  published  values  of  m^,  Mg  and  Mq  (seismic 
moment)  to  determine  spectral  scaling  relations  for  plate-margin  and 
mid-plate  earthquakes.  Following  Aki  (1967),  he  assumed  that  the  1-sec 
and  20-sec  period  amplitudes  of  the  far-fleld  spectrum,  after  correction 
for  attenuation,  are  proportional  to  m^  and  Mg,  respectively,  and  that 


14 


the  spectral  amplitude  at  very  long  periods  is  proportional  to  Mq. 
Nuttli  (1983a)  found,  by  trial  and  error,  sets  of  spectra  which,  on 
average,  satisfied  the  m^,  Mg,  Mq  combinations  for  mid-plate  and  plate- 
margin  earthquakes.  A  spectrum  of  the  mid-plate  spectra  consisted  sim¬ 
ply  of  a  linear  long-period  segment  of  constant  amplitude  (slope  zero) 
and  a  linear  short-period  portion  of  slope  two.  The  periods  at  which 

the  two  linear  segments  of  the  spectrum  intersect,  called  the  corner 

-4 

period  (T  ),  satisfied  the  relation  M  T  =  constant.  The  derived 
c  o  c 

plate-margin  spectra  were  more  complex,  because  for  seismic  moments 

23 

greater  than  about  10  dyne-cm  they  had  an  intermediate-period  linear 
segment  of  slope  one.  For  both  sets  of  spectra  there  is  the  same  one- 
to-one  correspondence  between  M_  and  M  values.  For  example,  for  Mc  = 
6,  the  spectral  amplitude  is  the  same  for  both  plate-margin  and  mid¬ 
plate  earthquakes.  However,  for  m^  (as  obtained  from  teleseismic  P-wave 
amplitudes),  a  value  of  0.4  (on  average)  had  to  be  added  to  the  of 
plate-margin  earthquakes  tc  result  in  the  same  spectral  amplitude  as 
that  corresponding  to  mid-plate  earthquakes.  That  is, 

mjj  (plate-margin  earthquakes)  +  0.4  =  m^  (mid-plate  earthquakes) 

This  equation  suggests  that  teleseismic  P-wave  amplitudes  of  plate- 
margin  earthquakes  suffer  a  2.5  times  (log-1 0.4)  larger  asthenosphere 
attenuation  in  the  source  region  than  mid-plate  earthquakes,  on  average. 
This  is  an  example  of  the  so-called  magnitude  bias  of  1-sec  period  P 
waves,  which  must  be  taken  into  consideration  when  m^  (teleseismic  P)  is 
used  to  estimate  explosion  yield. 

Assuming  values  of  average  fault  displacement  and  fault  rupture 
area,  Nuttli  (1983a)  showed  how  the  fault  rupture  length  and  average 
stress  drop  varied  with  magnitude  and  Mg,  and  seismic  moment  Mq,  for 
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average  mid-plate  and  plate-margin  earthquakes.  In  general,  for  a  given 
Mg  value,  plate-margin  earthquakes  have  lower  average  stress  drop  and 
greater  fault  rupture  length  than  mid-plate  earthquakes.  This  observa¬ 
tion  has  some  relevance  to  nuclear  explosion  detection.  For  mid-plate 
regions,  because  the  fault  rupture  lengths  of  moderate- sized  to  large 
earthquakes  are  relatively  small,  rather  insignificant  looking  geologi¬ 
cal  structures  can  produce  such  earthquakes.  Thus  for  mid-plate 
environments  it  is  more  difficult  to  identify,  from  surface  geological 
or  satellite  studies,  potential  earthquake  source  zones.  This  problem 
is  well  known  for  the  eastern  United  States,  and  often  produces  contro¬ 
versy  in  the  siting  of  nuclear  power  plants. 

Nuttli  (1983b)  added  more  data  and  did  more  extensive  study  of 
mid-plate  scaling  relations,  in  a  follow-up  of  the  earlier  paper 
(Nuttli,  1983a).  In  the  later  paper  he  concluded  that  the  spectra  of 
mid-plate  earthquakes  of  moment  greater  than  10  J  dyne-cm  required  an 
Intermediate- period  portion  of  unit  slope.  From  the  two  corner  periods 
the  fault  rupture  length  and  width  could  be  determined  as  a  function  of 
seismic  moment.  Also,  the  average  fault  displacement  could  be  calcu¬ 
lated,  using  equations  developed  by  Savage  (1972). 

Figure  1  compares  the  m^  versus  Mg  relation  for  mid-plate  and 
plate-margin  earthquakes,  as  determined  by  Nuttli  (1983a,  b).  The 
curves,  particularly  that  for  plate-margin  earthquakes,  indicate  a  more 
complicated  relation  than  the  generally  assumed  linear  one.  Also 
included  in  the  figure  is  the  explosion-earthquake  discriminant  curve  of 
Sykes  and  Evernden  (1982).  From  it  we  can  see  that  the  m^sMg  discrim¬ 
inant  is  more  effective  for  large  events  than  for  moderate-sized  ones 


5 

mb  (teleseismic  P) 


Fig.  1,  Relation  for  Mid-Plate  Earthquakes,  Plate-Margin  Earthquakes  and 

Explosions.  Included  are  the  data  for  two  New  Brunswick  earthquakes. 
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(m^  of  4.5  to  5.0) i  particularly  for  mid-plate  earthquakes.  Also  plot¬ 
ted  in  the  figure  are  the  m^Mg  values  (USGS,  1982)  of  the  9  January 
1982  New  Brunswick  earthquake  and  its  aftershock.  Both  practically  lie 
or  the  discriminant  curve  of  Sykes  and  Evernden. 

Chung  and  Bernreuter  (1981)  found  for  the  western  United  States 

that 

(teleseismic  P)  =  0.99  -  0.39 

Also  for  the  western  United  States,  Herrmann  and  Nuttli  (1982)  found 
that 

VLg)  .  ^ 

In  an  earlier  publication,  Nuttli  (1973)  defined  the  m.(Lg)  scale 
for  eastern  North  American  earthquakes  in  such  a  way  that,  for  that 
region, 

m^  (teleseismic  P)  =  m^Lg) 

From  these  three  relations,  and  if  we  assume  that  (teleseismic  P)  is 
reduced  by  0.4  units  because  of  abnormal  asthenosphere  attenuation,  we 
obtain  that  m^(Lg)  is  the  same  for  eastern  and  western  earthquakes  of 
the  same  1-sec  period  spectral  amplitude.  That  is,  ^(Lg)  is  an  esti¬ 
mate  of  the  1-sec  period  source  spectral  amplitude  that  is  Independent 
of  region.  The  conclusion  is  reasonable,  for  Lg  amplitudes  are  affected 
only  by  transmission  through  the  Earth's  crust.  The  conclusion  follows 
that  explosion  yield  estimates  based  on  m^(Lg)  values  should  exhibit  or 
have  no  magnitude  bias.  It  follows  that,  for  a  given  source  rock  type, 
suoh  as  hard  rock,  granite,  alluvium,  etc.,  there  will  be  an  m^tLg) 
versus  ylsld  ourve  that  will  be  applicable  to  all  continental  regions. 
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Preliminary  estimates  of  the  magnitude  bias  at  the  Nevada  (NTS), 
Bast  Kazakh  (USSR)  and  Sahara  (French)  test  sites,  relative  to  the 
eastern  United  States,  were  obtained  (Nuttli,  1982a,  b).  For  the  East 
Kazakh  site,  it  is  approximately  0.0  magnitude  units,  and  for  the  NTS 
and  Sahara  sites  it  is  0.4  units.  That  is,  if  eastern  North  America  is 
defined  to  have  zero  bias,  then  East  Kazakh  also  has  zero  bias,  whereas 
NTS  and  the  Sahara  sites  have  a  bias  of  0.4  units.  The  surprising 
result  is  the  large  bias  of  the  Sahara  test  site,  similar  to  that  of  the 
western  United  States  rather  than  to  a  shield  area. 
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Attenuation  of  Short  Period  Seismic  Crustal  Waves  in  Eurasia 

Otto  W.  Nuttli 

If  regional  seismic  waves  are  to  be  used  to  detect,  identify  and 
estimate  yields  of  explosions  in  the  USSR,  it  is  essential  that  their 
attenuation  be  accurately  known.  This  is  not  a  simple  task,  because  the 
attenuation  coefficients  vary  both  geographically  and  with  wave  fre¬ 
quency.  Under  Contract  F49620-C-79-0025  particular  attention  was  paid 
to  1-sec  period  Lg-wave  attenuation.  Values  of  this  parameter  were 
determined  for  paths  from  the  East  Kazakh  test  site  to  the  WWSSN  seismo¬ 
graph  stations  in  Scandinavia  and  southern  Asia,  for  paths  across  Iran 
and  for  paths  from  other  explosions  and  from  earthquakes  in  the  USSR  to 
Scandinavian  and  southern  Asian  stations.  In  addition,  the  attenuation 
coefficients  of  1-sec  period  Pn,  Pg  and  Sn  waves,  as  well  as  3-sec 
period  Lg  waves,  were  determined  for  paths  across  Iran. 

Several  methods  were  used  to  determine  the  values  of  the  attenua¬ 
tion  coefficients  of  the  crustal  waves.  When  sufficient  data  were 
available,  time-domain  amplitudes  were  plotted  against  epicentral  dis¬ 
tance  and  fitted  by  theoretical  curves  whose  shapes  are  a  function  of 
the  attenuation  coefficient  value.  This  was  the  method  employed  by 
Nuttli  (1980a)  for  a  study  of  1-sec  and  3-sec  period  Lg  waves  and  of  1- 
sec  period  Pn,  Pg  and  Sn  waves  across  Iran.  This  method  gives  average 
values  over  broad  areas  or  over  long  profiles.  A  second  method 
developed  by  Herrmann  (1980),  in  an  adaptation  of  Aki's  (1969)  theory 
that  modeled  the  Lg  coda  as  scattered  waves,  was  used  to  determine  the 
average  attenuation  coefficient  for  Lg  waves  of  a  given  period  for  indi¬ 
vidual  paths  from  source  to  receiver.  This  method  has  the  advantage  of 
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giving  the  attenuation  for  a  particular  path,  which  is  necessary  when 
the  attenuation  varies  rapidly  over  an  area.  It  is  especially  useful 
when  there  are  many  seismic  events  that  are  located  in  a  small  source 
region,  such  as  a  test  site,  and  are  recorded  at  the  same  group  of  sta¬ 
tions.  A  second  advantage  of  the  Lg-coda  method  is  that  the  attenuation 
can  be  determined  from  the  variation  of  the  wave  period  in  the  Lg  coda 
with  travel  time,  and  thus  is  independent  of  amplitude  measurements. 
There  are,  however,  also  some  limitations  or  problems  with  the  coda 
method  for  determining  attenuation.  The  principal  one  is  that  the  value 
of  the  attenuation  coefficient  is  frequency  dependent,  and  the  coda  data 
often  are  not  available  over  a  sufficient  range  of  periods  to  determine 
this  frequency  dependence  precisely.  Other  methods  then  have  to  be 
resorted  to  in  order  to  fix  the  value  of  in  the  relation 

Q(f)  +  Q0f^  (Mitchell, 
where  f  is  wave  frequency,  ^  is  taken  to  be  constant  over  the  range  of 
about  0.2  to  10  sec,  and  Q  is  the  specific  quality  factor,  related  to 
the  anelastic  attenuation  coefficient,  y,  by  the  equation 

y(f)  +  *f/UQ(f) 

where  0  is  the  group  velocity. 

For  Iran,  Nuttli  (1980a)  found  the  value  of  the  anelastic  attenua¬ 
tion  coefficient  of  1-sec  period  Pg,  Sn  and  Lg  waves  to  be  0.0045  km"1, 
approximately  the  same  value  that  was  found  for  Lg  in  California  by 
Herrmann  (1980).  For  3-sec  period  Lg  waves  propagating  across  Iran,  the 
value  of  >  was  found  to  be  0.003  km"1.  The  corresponding  approximate 
relation  for  Lg  waves  in  Iran  is 


Q(f)  =  200f0*6 

Making  use  of  the  fall-off  of  amplitude  with  epioentral  distance  in 
Iran,  Nuttli  (1980a)  developed  formulas  for  m^,  valid  for  distances  of 
400  to  1500  km,  that  utilize  the  amplitude  of  1-sec  period  Pg  or  Lg 
waves  recorded  by  vertical-component  seismographs.  The  m^  obtained  from 
these  formulas  was  defined  to  be  the  same  as  that  obtained  from  telese- 
ismic  P-wave  amplitudes  for  Iranian  earthquakes. 

Also  it  was  found  that  t.»e  excitation  of  1-sec  period  Lg  waves  for 
an  m^  (teleseismic  P)  5.0  earthquake  in  Iran  is  the  same  as  in  Califor¬ 
nia,  and  approximately  0.4  magnitude  units  greater  than  for  a  similar- 
sized  earthquake  in  eastern  North  America.  An  alternative  way  of  stat¬ 
ing  this  is  that  the  teleseismic  P-wave  amplitudes  of  seismic  events 
occurring  either  in  Iran  or  in  California  suffer  an  anomalous  0.4  magni¬ 
tude  unit  loss  of  amplitude,  relative  to  eastern  North  America,  in  pro¬ 
pagating  through  the  asthenosphere  beneath  the  source.  Thus,  the  so- 
called  magnitude  bias  of  Iran.. in  and  Californian  events,  relative  to 
eastern  North  America,  is  0.4  units. 

Nuttli  (1981)  gave  values  of  1-sec  period  Lg  attenuation  coeffi¬ 
cients  for  paths  from  16  explosions  in  western  Russia,  the  Urals, 
western  Kazakhstan,  and  western  Siberia  and  from  6  earthquakes  in 
eastern  Caucasus  to  13  WWSSN  stations  in  southern  Asia,  Scandinavia  and 
Germany.  These  values  were  obtained  by  application  of  Herrmann's  (1980) 
coda  method. 

The  coda  method  also  was  used  by  Nuttli  (1982)  to  determine  Qq  and 
$  values  for  paths  from  the  East  Kazakh  test  site  to  southern  Asian  and 
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Scandinavian  stations.  The  values  found  were 


Station 

Country 

Qo 

§ 

KBL 

Afghanistan 

360 

0.6 

KEV 

Finland 

580 

0.4 

KON 

Norway 

700 

0.4 

MKI 

Iran 

380 

0.5 

NDI 

India 

300 

0.6 

NIL 

Pakistan 

380 

0.6 

NUR 

Finland 

580 

0.4 

QUE 

Pakistan 

280 

0.6 

SHL 

India 

340 

0.6 

0ME 

Sweden 

620 

0.4 

In  an  earlier  report  (Nuttli,  1980b),  Qq  values  were  estimated  for 
the  Lg  waves  from  a  number  of  explosions  and  earthquakes  throughout 
western  and  central  Eurasia  recorded  at  the  same  stations  listed  above. 
In  this  early  study,  the  value  of  (j  was  assumed  to  be  0.2  for  all  the 
source- to-station  paths.  As  a  consequence,  because  this  value  of  (j  is 
too  small,  all  the  Qq  values  given  in  Nuttli  (1980b)  were  underes¬ 
timated,  more  so  for  the  paths  to  the  southern  Aslan  stations  than  to 
those  in  Scandinavia.  The  more  recently  obtained  ^  values  of  0.4  to  0.6 
(see  table  above)  are  consistent  with  (j  values  of  0.5  for  S  waves  in 
central  Asia  observed  by  Rautian  and  Khlaturin  (Aki,  1980),  of  0.6  for  S 
waves  in  the  Kurils  observed  by  Fedotov  and  Boldyrev  (Aki,  1980)  and  of 
0.6  and  0.8  for  S  waves  observed  by  Aki  (1980)  in  Japan. 
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EXPLOSION  YIELD  ESTIMATES 
USING  REGIONAL  Lg  WAVES 
Otto  W.  Nuttli 

The  conventional  seismic  methods  for  estimating  the  yields  of 
underground  nuclear  explosions  make  use  of  body-wave  magnitude,  m^,  and 
surface-wave  magnitude,  Mg.  This  is  done  by  .setting  up  a  correlation 
between  each  of  the  magnitudes  and  known  explosion  yields. 

Both  methods  suffer  from  certain  complicating  factors.  For  exam¬ 
ple,  values  are  affected  by  the  type  of  rock  in  the  source  region. 

In  addition,  a  so-called  magnitude  bias  is  caused  by  anelastic  attenua¬ 
tion  of  1-Hz  P-wave  amplitudes  in  the  asthenosphere.  The  value  of  this 
bias  term  varies  from  one  geographic  region  to  another.  Also,  focusing 
and  defocusing  of  short-period  P  waves  by  lateral  crustal  variations  at 
a  test  site  can  produce  anomalous  m^  values.  Empirical  correlation 
between  and  yield  should  be  corrected  for  these  factors,  which  may 

not  be  well  known  for  different  test  sites. 

Because  the  waveguide  for  20-sec  period  surface  waves  is  rather 
uniform  over  much  of  the  Earth,  the  problems  that  are  encountered  for  P 
waves  are  not  so  serious  for  surface  waves.  But  there  are  two  other 
types  of  problems.  The  first  is  that  20-sec  period  waves  do  not  have 
sufficient  amplitude  to  be  used  for  estimating  yields  of  less  than  10 
kt.  When  large  amplitude  surface  waves  of  period  about  10  sec  are  used 
to  overcome  this  limitation,  it  is  found  that  they  exhibit  large  lateral 
variations.  The  second,  and  more  serious  problem,  is  that  nuclear 
explosions  can  release  tectonic  strain  and  produce  an  earthquake- type 
radiation  pattern  superposed  on  the  isotropic  radiation  field  of  an 
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explosion.  Depending  upon  the  distribution  of  seismograph  stations,  the 
net  result  oan  be  either  to  increase  or  decrease  the  estimate  of  the 
explosion  yield. 

As  a  result  of  these  problems,  and  Mg  values  can  sometimes  pro¬ 
duce  significantly  different  values  of  the  estimated  yield  of  a  particu¬ 
lar  explosion. 

A  third  type  of  wave,  the  higher-mode  surface  wave  Lg,  shows  prom¬ 
ise  of  providing  an  independent  estimate  of  yield.  It  has  certain 
advantages  and  disadvantages  with  respect  to  methods  using  P  and  surface 
waves.  An  obvious  disadvantage  is  that  it  only  can  be  applied  for  con¬ 
tinental  paths  from  source  to  receiver,  because  1-sec  period  Lg  waves  do 
not  propagate  across  oceanic  paths.  A  second  disadvantage  i3  that  it  is 
limited  to  regional  distances,  which  in  general  restricts  its  applica¬ 
tion  to  distances  of  less  than  1000  tan  for  neotectonic  areas  and  about 
5000  tan  for  ancient  shield  areas. 

The  advantages  of  using  Lg  amplitudes  for  yield  estimates  are  that 
the  method  can  be  applied  to  small  events  (less  than  1  kt),  that  there 
appears  to  be  no  magnitude  bias,  such  as  observed  for  P  waves,  and  that 
apparently  the  method  is  not  affected  by  the  earthquake- type  radiation 
produced  by  the  release  of  tectonic  strain,  so  that  a  single  calibration 
curve  can  be  used  for  estimating  yield  for  all  geographic  regions.  How¬ 
ever,  this  latter  point  needs  to  be  tested  by  further  study  and  docu¬ 
mented  before  it  can  be  accepted. 

For  the  past  ten  years  Lg  amplitudes  have  been  used  routinely  in 
eastern  North  America  for  estimating  body-wave  magnitude,  called  a*^(Lg) 


(Nuttli,  1973).  In  general,  1-sec  Lg  amplitudes  for  eastern  earthquakes 
give  m^  values  equal  to  m^  from  1-sec  period  P-wave  amplitudes,  as 
Intended.  They  scale  In  the  same  way,  because  both  are  a  measure  of  the 
1-sec  period  spectral  amplitude.  The  Lg  method  has  proven  useful 
because  it  can  be  applied  to  events  of  m^  as  small  as  2.5,  when  WWSSN 
seismograms  are  used. 

Fig.  1  illustrates  how  Lg  amplitudes  sire  used  to  determine  magni¬ 
tude.  The  maximum  sustained  amplitude  at  a  period  near  1  sec  and  a 
group  velocity  near  3.5  km/ sec  is  measured  on  the  vertical-component 
seismograms,  and  converted  to  ground  motion.  The  amplitude  then  is 
extrapolated  back  to  10-km  epicentral  distance  according  to  the  formula 
(adapted  from  Ewing  at  si. ,  1957) 

A(10  km)  =  A(£)  (/V10)1/3  [sin(^yi11)/sin(10/111)]1/2  exp[X(/\-  10)] 
where  is  epicentral  distance  in  kilometers.  The  quantity  y  is  called 
the  coefficient  of  anelastic  attenuation.  It  is  given  by 

X(f)  =  ff/UQ( f ) 

where  f  is  wave  frequency,  U  is  the  group  velocity  of  the  wave  whose 
amplitude  was  read,  and  Q  is  an  apparent  specific  quality  factor.  Q  is 
assumed  to  be  frequency  dependent  according  to  the  relation  proposed  by 
Mitchell  (1980) 

Q(f)  =  Q0f^ 

where  Qq  is  the  value  of  Q  at  1-sec  period. 

For  earthquakes  in  eastern  North  America,  an  (tel.  P)  =  5.0 

corresponds  to  an  A(10)  value  of  the  vertical  component  of  Lg  of  110 
microns.  Thus,  on  the  condition  that  m^(Lg)  and  m^  (tel.  P)  scale  in 
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y(f)  ■  *  f  /U(f)  0(f)  and  Q(f)  ?  Q#  f 5 

FOR  m^«5.0t  A (I0)z *  110 microns 
«»b(Lg)»  3.0  |og  [A(I0)/|I0] 


Fig.  1.  Methodology  for  Using  Lg  Amplitudes  and  Frequencies  to 
Determine  mj>(Lg) 


m^(Lg)  =  5.0  +  log[A( 10)/1 10] 

where  A(10)  is  the  extrapolated  value  of  the  1-sec  period,  vertical- 


component  Lg  ground  motion  (in  microns)  at  a  distance  of  10  km. 

In  the  early  work  X  (the  coefficient  of  anelastic  attenuation)  was 
obtained  by  measuring  the  decrease  of  amplitude  with  epicentral  dis¬ 
tance.  This  method  was  used  to  obtain  Y  values  for  1-sec  period  Lg 
waves  in  eastern  North  America  and  in  Iran,  and  of  shorter  period  waves 
in  the  Mississippi  Valley  and  in  California.  Later  we  observed  that  the 
coda  theory  of  Aki  and  Chouet  (1975),  as  adapted  by  Herrmann  (1980), 
gave  the  same  values  of  X  as  the  fall-off  of  amplitude  with  distance. 
The  coda  method  is  ideally  suited  for  nuclear  explosion  studies,  for  it 
gives  the  average  Qq  and  ^  values  for  the  path  from  source  to  receiver. 
Fig.  2  shows  an  example  of  the  fit  of  some  theoretical  curves  to  the 
data  for  paths  from  NTS  to  the  stations  DOG  (  A  Z  450  km) ,  BKS  (  ;  550 
km) ,  and  TOC  (  Z  700  km) . 

Measurements  of  coda  waves  for  sources  in  the  various  parts  of  the 
NTS  (e.g.  Yucca  Flats  in  the  east  central  and  Pahute  Mesa  in  the 
northwest  portions  of  the  Test  Site)  suggest  that  the  Qq  and  (}  values 
for  the  paths  to  the  stations  DOG,  BKS  and  TOC  do  not  vary  with  source 
location.  Thus,  lateral  variations  in  crustal  structure  over  the  NTS  do 
not  appear  to  affect  1-sec  Lg  amplitudes,  in  contrast  to  what  is 
observed  for  1-sec  P-wave  amplitudes. 

Springer  and  Kinnaman  (1971,  1975)  published  yields  for  a  number  of 
explosions  at  NTS,  along  with  information  about  location,  origin  time, 
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Fig.  2.  Qc  and  S  Value*  for  Paths  from  Yucca  Flats, NTS  to  BKS,  DUG  and  TUC 
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source  rock  type  and  depth  of  the  shot.  By  examination  of  the  Lg  waves 
from  these  explosions  we  found  that  DOG,  BKS  and  TUC  were  the  best  WWSSN 
stations  available  for  this  study.  However,  they  were  not  ideal.  For 
example,  TOC  operates  at  a  magnification  of  200,000  at  1-sec  period, 
which  results  in  the  Lg  waves  going  off  the  edges  of  the  seismogram  for 
explosions  of  about  50  kt  and  greater.  DOG,  which  is  closer  to  NTS  than 
TOC,  normally  operates  at  a  1-sec  magnification  of  200,000  or  400,000. 
Therefore,  the  motion  "whited  out"  or  went  offscale  for  events  of  yield 
as  low  as  20  kt,  except  when  the  gain  was  reduced  before  some  of  the 
larger  explosions. 

Fig.  3  shows  the  ff^(Lg)  versus  log  yield  data  and  curve  for  explo¬ 
sions  of  announced  yield  in  hard  rock  at  NTS.  The  curve  is  a  quadratic 
least-squares  fit  to  the  data,  given  by  the  equation 

mb  =  3-901  +  1.185  log  Y  -  0.102  (log  Y)2 
where  Y  is  the  announced  yield  in  kilotons.  Considering  log  Y  to  be 
error-free,  the  standard  deviation  of  is  0.077.  The  slope  of  the 
curve  at  small  yields  (10  kt)  is  0.98,  very  close  to  unity  as  is  to  be 
expected  from  the  spectral  scaling  relations  for  explosions  given  by 
Mueller  and  Murphy  (1971).  At  1000  kt  the  slope  has  a  value  of  0.57, 
which  also  is  in  agreement  with  the  spectra  scaling  relations. 

If  SHOAL  and  PILE  DRIVER,  which  were  shot  in  granite,  are  not 
inoluded  in  the  quadratic  least-squares  fit,  the  resulting  equation  is 

=  3.901  +  1.155  log  Y  -  0.091  (log  Y)2 
with  a  standard  deviation  of  m^  of  0.066.  The  resulting  curve  is  so 
similar  to  that  of  Fig.  3  that  it  is  not  replotted.  The  slope  at  a 
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yield  of  10  kt  is  0.97,  and  at  1000  kt  is  0.61. 

Pig.  4  gives  the  m^(Lg)  versus  log  yield  data  and  curve  for  explo¬ 
sions  in  alluvium  at  NTS.  Although  a  quadratic  equation  was  used  to  fit 
the  data,  the  curve  almost  is  linear.  Its  equation  is 

=  4.041  ♦  0.726  log  Y  -  0.009  (log  Y)2 
with  a  standard  deviation  of  0.048  for  m^,  assuming  log  Y  is  error-free. 
For  Y  s  10  kt,  m^Lg)  For  a  source  in  alluvium  is  approximately  0.2 
units  less  than  for  a  -source  in  hard  rock.  At  Y  =  1000  kt,  the  differ¬ 
ence  in  m^  values  increases  to  about  0.4  units. 

Table  1  shows  the  differences  between  the  m^  values  obtained  from 
teleseismlc  P-wave  amplitudes  and  regional  Lg-wave  amplitudes.  On  aver¬ 
age,  m^(Lg)  at  NTS  is  0.31  units  larger  than  ^(P).  From  earthquake 
magnitude  studies,  Chung  and  Bernreuter  (1981)  and  Herrmann  and  Nuttli 
(1982)  concluded  that  m^P)  is  reduced  by  about  0.4  units  because  of 
anelastic  attenuation  of  P  waves  in  the  asthenosphere  beneath  the 
western  United  States.  Therefore,  for  earthquakes,  m,o(Lg)  should  be  0.4 
units  larger  than  m^P).  The  difference  between  0.4  and  0.31  may  be 
considered  a  measure  of  the  discrimination  potential  of  1-sec  period  Lg 
waves,  i.e,,  how  much  the  amplitude  of  explosion-generated,  1-sec  period 
Lg  waves  differs  from  that  of  earthquake-generated  Lg  waves. 
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TABLE  1 


Date 

Name 

NTS  HARE 

Rock  Type 

i  ROCK 

^(Lg) 

SITES 

Station 

(P-vfec) 

Announced 

(kt) 

10/05/62 

Mississippi 

tuff-above 

5.78 

BKS 

__ 

110 

09/13/63 

Bilby 

tuff-below 

6.06 

BKS 

— 

235 

10/26/63 

Shoal 

granite 

5.19 

BKS, TUC 

— 

12 

11/05/64 

Handcar 

dolomite 

5.04 

BKS, DUG, TUC 

4.8 

12 

12/16/64 

Mudpack 

tuff-above 

4.32 

BKS, TUC 

2.7 

04/14/65 

Palanquin 

rhyolite 

4.67 

BKS,DUG,TUC 

— 

4.3 

02/24/66 

Rex 

tuff-below 

5.28 

BKS, TUC 

5.0 

19 

04/14/66 

Duryea 

rhyolite 

5.66 

BKS.DUG.TUC 

5.4 

65 

05/06/66 

Chartreuse 

rhyolite 

5.76 

BKS, DUG 

5.4 

70 

05/27/66  Discus  Thrower  tuff-above 

5.25 

BKS, TUC 

5.0 

21 

06/02/66 

Pile  Driver 

granite 

5.82 

BKS 

5.6 

56 

06/30/66 

Half  Beak 

rhyolite 

6.25 

BKS, TUC 

6.1 

300 

12/20/66 

Greeley 

tuff-below 

6.49 

BKS, DUG 

6.3 

825 

05/20/67 

Commodore 

tuff-below 

6.00 

BKS 

5.8 

250 

05/23/67 

Scotch 

tuff-below 

5.96 

BKS 

5.7 

150 

05/26/67 

Knickerbocker 

rhyolite 

5.75 

BKS, DUG 

5.4 

71 

01/26/68 

Cabriolet 

rhyolite 

4.34 

DUG, TUC 

— 

2.3 

03/12/68 

Buggy  I 

basalt 

4.63 

BKS, DUG, TUC 

— 

5.4 

04/26/68 

Boxcar 

rhyolite 

6.58 

BKS 

6.2 

1300 

12/08/68 

Schooner 

tuff-above 

5.32 

BKS, DUG, TUC 

4.8 

30 

12/19/68 

Benham 

tuff-below 

6.61 

BKS, DUG 

6.3 

1100 

10/29/69 

Calabash 

tuff 

5.88 

BKS, DUG 

5.6 

110 

03/26/70 

Handley 

tuff-below 

6.58 

BKS, DUG 

6.4 

>1000 

05/26/70 

Flask 

tuff-above 

5.86 

BKS, DUG 

5.5 

105 

10/17/70 

Carpetbag 

tuff-below 

6.14 

BKS, DUG 

5.8 

220 

07/08/71 

Miniata 

tuff-above 

5.94 

BKS, DUG, TUC 

5.5 

80 

04/26/87 

Starwort 

tuff-below 

5.73 

BKS, DUG, TUC 

5.6 

85 

10/09/64 

Pan 

NTS  ALLUVIUM  SITES 

alluvium  5.23  BKS, DUG, TUC 

4.8 

38 

12/16/64 

Parrot 

alluvium 

4.06 

DUG, TUC 

— - 

1 .2 

02/16/65 

Merlin 

alluvium 

4.72 

BKS, DUG, TUC 

— 

10 

06/11/65 

Pet:  1 

alluvium 

4.19 

DUG, TUC 

— 

1.2 

05/05/66 

Cyclmen 

alluvium 

4.86 

NKS,DUG,TUC 

13 

09/21/67 

Marvel 

alluvium 

4.30 

DUG, TUC 

— 

2.2 

05/14/68 

Pommord 

alluvium 

4.09 

DUG , TUC 

— - 

1.4 

09/26/72 

Delphinium 

alluvium 

4.87 

BKS, TUC 

4.4 

15 

tBb(lg)  -  Wage  =  °’31  *  °*10  (n  =  23> 
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Regional  Lg  waves  from  East  Kazakh  (Degelen  Mountain  and  Shagan 
River)  explosions  are  observed  at  a  number  of  WWSSN  stations  in  Scandi¬ 
navia  and  southern  Asia.  They  can  be  used  to  determine  m^Lg).  Assum¬ 
ing  that  the  relation  between  m.  (Lg)  and  yield  is  the  same  for  hard-rock 
explosions  at  NTS  and  East  Kazakh,  the  yields  of  East  Kazakh  explosions 
can  be  estimated. 

Table  2  shows  the  Q  values,  the  upper  and  lower  limits  of  Q 
obtained  by  the  coda  data,  and  the  change  in  Qq  required  to  cause  a  O.i 
unit  change  in  m.(Lg),  or  a  25$  change  in  estimated  yield  for  the  small 
yields  (<100  kl)  and  a  50$  change  in  estimated  yield  for  the  large 
yields  (>400  kt). 

Table  3  gives  preliminary  estimates  of  m.  (Lg)  and  yield  for 
selected  East  Kazakh  events.  Of  the  42  explosions  listed  in  Table  3, 
the  average  Bi^(Lg)  -  m.  (tel.  P)  =  -0.14.  For  NTS,  the  average  differ¬ 
ence  was  +0.31  units.  The  difference  between  -0.14  and  +0.31 ,  or  0  “j 
units,  could  be  caused  by  a  different  relation  between  m^(Lg)  and  yield 
for  NTS  and  East  Kazakh  explosions  (different  excitations  of  Lg)  or  by 
an  anomalous  absorption  of  P  amplitudes  at  NTS  (relative  to  East 
Kazakh) .  The  latter  explanation  appears  more  likely,  as  it  would 
approximately  agree  with  Chung  and  Bernreuter's  (1981)  and  Herrmann  and 
Nut til's  (1982)  estimate  of  anomalous  P-wave  attenuation  in  the  astheno- 
sphere  beneath  the  western  United  States. 
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TABLE  2 

SOME  PARAMETERS  FOR  EAST  KAZAKH  SOURCES 


Station 

Average 

Distance  (ten) 

*o 

$ 

£flo» 

KBL 

1900 

360  (340-380) 

0.6 

17 

KEY 

3500 

580  (550-620) 

0.4 

28 

KON 

4300 

700  (650-760) 

0.4 

42 

MHI 

2100 

380  (350-410) 

0.5 

18 

NDI 

2400 

300  (270-320) 

0.6 

13 

NIL 

1900 

380  (360-410) 

0.6 

21 

NUR 

3500 

580  (520-650) 

0.4 

31 

QUE 

2200 

280  (260-310) 

0.6 

10 

SHL 

2900 

340  (320-350) 

0.6 

11 

UME 

3700 

620  (610-640) 

0.4 

30 

s  change  In  Qq  to  causa  0-10  unit  change  in  Df,(Lg). 


j 
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TABLE  3 

m^(Lg)  AND  PRELIMINARY  ESTIMATED  YIELDS  FOR 
SELECTED  EAST  KAZAKH  EXPLOSIONS 


Date 

Location 

®lj(Lg) 

Number  of 
Stations 

mb(P)-ISC 

Estimated 
Yield  (kt) 

01-15-65 

Shagan 

6.02 

2 

5.8 

166 

11-02-72 

S 

6.04 

1 

6.1 

178 

12-10-72 

S 

6.06 

2 

6.0 

186 

04-19-73 

Degelen 

5.18 

2 

5.4 

17 

07-10-73 

D 

5.23 

2 

5.2 

19 

07-23-73 

S 

6.22 

6 

6.1 

309 

10-26-73 

D 

4.99 

1 

5.2 

10 

12-14-73 

S 

5.77 

3 

5.8 

79 

01-30-74 

D 

5.27 

2 

5.4 

21 

05-16-74 

D 

4.94 

2 

5.2 

9.3 

05-31-74 

S 

5.75 

4 

5.9 

74 

07-10-74 

D 

4.90 

2 

5.2 

8.5 

09-13-74 

D 

4.97 

1 

5.2 

10 

10-16-74 

S 

5.23 

3 

5.5 

19 

12-16-74 

D 

4.20 

2 

4.8 

1.8 

12-27-74 

S 

5.68 

6 

5.6 

62 

02-20-75 

D 

5.41 

3 

5.7 

30 

03-11-75 

D 

5.18 

2 

5.4 

17 

04-27-75 

S 

5.38 

4 

5.6 

28 

06-08-75 

D 

5.11 

2 

5.5 

14 

08-07-75 

D 

4.96 

2 

5.2 

9.8 

10-29-75 

S 

5.45 

2 

5.8 

33 

12-25-75 

S 

5.82 

4 

5.7 

91 

01-15-76 

D 

5.01 

2 

5.2 

11 

06-09-76 

S 

5.26 

3 

5.3 

20 

07-04-76 

S 

5.75 

2 

5.8 

74 

08-26-76 

S 

5.54 

5 

5.8 

42 

12-07-76 

S 

5.64 

2 

5.9 

55 

12-30-76 

D 

5.02 

1 

11 

09-05-77 

S 

5.46 

2 

5.8 

34 

11-30-77 

S 

5.63 

6 

6.0 

54 

07-05-78 

S 

5.61 

1 

5.8 

51 

08-29-78 

S 

5.70 

3 

5.9 

65 

09-15-78 

/  S 

5.78 

4 

6.0 

81 

11-04-78 

s 

5.59 

2 

5.6 

48 

11-29-78 

s 

5.70 

4 

5.3 

65 

06-23-79 

s 

5.98 

2 

6.2 

148 

07-07-79 

s 

5.88 

2 

5.8 

110 

08-04-79 

s 

5.97 

4 

6.1 

141 

08-18-79 

s 

6.01 

5 

6.1 

162 

09-14-80 

s 

6.07 

5 

6.2 

195 

avg.  m^(Lg)  -  m^P)  =  0.15  i.  0.19  (N  =  41) 
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Multi-mode  Surface-wave  Propagation 

and  Qp  the  Continental  Crust 
by 

Brian  J.  Mitchell 

Recent  studies  have  indicated  that  Q  structure  varies  significantly 
from  region  to  region  throughout  continents  as  well  as  oceans.  Mitchell 

(1975)  showed  that  values  in  the  upper  crust  over  a  broad  region  of 
western  North  America  were  about  half  as  large  as  those  values  in  the 
upper  crust  of  eastern  North  America. 


In  order  to  understand  the  regional  variations  of  detail  it 
is  necessary  to  study  the  attenuation  of  seismic  waves  over  regions 
which  are  as  small  as  possible.  A  single  source-single  station  method 
has  been  developed  in  order  to  study  the  anelastic  properties  of  the 
crust  along  relatively  short  paths  (Cheng  and  Mitchell,  1981).  By  util¬ 
izing  short  paths,  effects  of  lateral  variation  of  Q  along  the  path,  as 
well  as  the  effects  of  refraction,  reflection,  and  multi- pathing  of  sur¬ 


face  waves  can  be  minimized.  We  have  now  applied  the  method  to  three 


regions  in  North  America,  the  eastern  United  States,  the  Colorado  Pla¬ 
teau,  and  the  Basin-and-Range  Province,  as  well  as  to  a  region  of  Iran 
and  Turkey  and.  the  Barents  shelf. 


Method 


Because  of  the  overlap  In  the  group  arrival  times  of  the  several 
higher  modes  which  comprise  the  phase  Lg  in  continental  regions,  it  is 
usually  not  possible,  using  single  station  recordings,  to  separate  those 
modes  for  individual  study.  To  avoid  that  problem,  Cheng  and  Mitchell 
(1981)  introduced  the  concept  of  multi-mode  spectra.  Using  that  con¬ 
cept,  theoretical  spectra,  consisting  of  the  contributions  of  several 
higher  modes,  can  be  compared  to  observed  spectra  corresponding  to  the 
composite  of  higher  modes  which  form  the  phase  Lg  on  long-period  seismo¬ 
grams.  An  example  appears  in  Figure  1 . 

There  are  four  assumptions  which  are  inherent  in  the  multi-mode 
method:  (1)  that  the  shape  of  the  source  spectrum  of  the  seismic  event 
can  be  estimated  within  certain  bounds,  (2)  that  velocity  structure 
along  the  path  of  propagation  is  approximately  known,  (3)  that  the 
mechanism  which  causes  seismic  waves  to  attenuate  affects  both  the  fun¬ 
damental  mode  and  higher  modes  in  the  same  way  and  that  the  attenuation 
can  be  described  by  exp[-  /Up^Qp]  where  Ur  is  the  group  velocity,  and  Qp 

is  the  Q  value  for  a  given  Rayleigh  mode  at  period  T,  and  (4)  that  Qn  is 

P  X 

independent  of  frequency. 

Assumption  (1)  requires,  in  the  case  of  earthquakes,  that  the  focal 
depth  and  fault-plane  solution  are  known,  and  since  we  usually  assume 
the  source  to  be  a  step- function  in  time,  that  the  observed  spectra  do 
not  depart  greatly  from  those  predicted  for  such  a  source.  For  explo¬ 
sion  sources,  we  should  take  account  of  a  small  maximum  in  the  source 
spectrum  at  short  periods  as  discussed  later.  The  absolute  level  of  the 
source  spectrum  need  not  be  known,  since  we  can  attempt  to  fit  the 


ler-mode  spectra  and  plus  signs  denote  "observed"  fundamental-mode 
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fundamental-  and  higher-mode  spectral  data  by  adjusting  seismic  moment, 
as  well  as  the  value  of  Q  in  the  upper  crust. 

Assumption  (2)  usually  poses  no  difficulties  since  we  can  use  the 
group  velocities  observed  for  the  analyzed  seismograms  to  determine  a 
crustal  velocity  model.  In  our  inversion,  we  neglect  the  dependence  of 
surface  wave  velocities  on  Q  structure  in  the  earth  and  invert  for  Q 
structure  alone,  assuming  our  velocity  model  is  reasonably  well  deter¬ 
mined. 


The  validity  of  assumption  (3)  is  the  most  difficult  of  our  assump¬ 
tions  to  assess.  It  is  possible,  for  instance,  that  the  energy  in  one 
mode  might  be  reflected  or  converted  to  other  types  of  waves  more  easily 
than  that  of  another  mode.  For  this  reason,  we  attempt  to  use  rela¬ 
tively  short  paths  across  regions  which  are  devoid  of  obvious  major 
lateral  changes.  In  addition,  we  attempt  to  use  several  paths  in  dif¬ 
ferent  locations  and  directions  and  to  find  a  model  which  explains  all 
of  the  data  as  well  as  possible. 

Assumption  (4)  may  not  be  correct  since  Qp  appears  to  increase  with 
decreasing  period,  at  least  at  periods  shorter  than  2  or  3  seconds 
(Mitchell,  1980).  As  shown  in  that  study,  however,  it  is  not  likely, 
using  surface  wave  data  and  presently  available  methods,  that  we  could 

detect  a  frequency  dependence  of  Qp  at  periods  between  about  3  and  50 
seconds.  In  the  absence  of  data  at  shorter  periods,  therefore,  we  must 
restrict  our  study  to  the  determination  of  frequency-independent  models. 

For  the  regions  discussed  below,  we  make  yet  another  assumption, 
that  the  value  of  Qp  the  lower  crust  is  much  higher  than  that  in  the 
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upper  crust.  This  assumption  is  based  both  upon  previous  determinations 

of>  in  the  central  United  States  where  the  best  short-period  surface 
wave  attenuation  data  have  been  obtained,  and  upon  the  idea  that  fluid- 

filled  cracks  and  pore-space  in  the  crust,  which  are  likely  contributors 
to  crustal  attenuation  close  at  mid-crustal  depths.  Although  we  have 


assumed  that  100O  in  the  lower  crust,  we  shall  see  that  our 

observed  spectra  are  not  very  sensitive  to  lower  crustal  values. 


The  multi-mode  method  considers  the  amplitude  spectrum  of  the  com¬ 
bined  higher  modes,  as  well  as  the  spectrum  of  the  fundamental  mode.  If 

Qp  in  the  lower  crust  is  much  higher  than  Qp  in  the  upper  crust,  then 
the  fundamental-mode  spectrum  is  more  strongly  affected  than  the 
higher-mode  spectrum  by  changes  of  Qp  in  the  upper  crust  (cheng  ^ 
Mitchell,  1981).  There  is,  of  course,  some  effect  of  the  shapes  of  both 

spectra  by  changes  in  Given  the  scatter  in  the  data,  however,  it  is 
not  possible  to  assign  significance  to  small  changes  in  spectral  shape. 
We  therefore  look  at  two  things:  the  ratio  of  the  fundamental-mode  spec¬ 
trum  to  that  of  the  higher  modes  and  the  gross  shape  of  the  spectra.  In 


situations,  such  as  those  discussed  below,  where 


values  in  the  lower 


crust  are  higher  than  those  in  the  upper  crust,  the  higher-mode  spectrum 


serves  mainly  as  a  reference  to  which  the  fundamental-mode  spectrum  is 


compared.  Of  course,  large  differences  between  the  shapes  of  the 


observed  and  theoretical  spectra  would  lead  us  to  try  other  Q  models. 


In  later  paragraphs  we  will  see  that  both  the  ratios  and  shapes  of  the 
spectra  are  primarily  sensitive  to  the  values  in  the  upper  crust. 


We  are  limited  in  our  interpretation  by  two  factors,  (1)  the  rela¬ 
tively  poor  resolution  obtainable  using  seismic  wave  amplitudes,  even  if 
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we  have  high-quality  data,  and  (2)  the  scatter  in  our  data  due  to  depar¬ 
tures  from  ideal  sources  and  ideal  propagation  conditions.  We  therefore 
attempt  to  find  only  very  simple  Q  models  of  the  crust.  These  models 
are  much  cruder  than  velocity  models  which  can  be  obtained  using  modern 
inversion  methods.  Our  purpose  is,  however,  to  look  only  for  major 

variations  of  Qp  t^e  upper  crust  of  various  regions  of  the  world. 
Fortunately,  the  regional  variation  of  Qp  at  those  depths  can  be  quite 
large  and  the  differences  from  region  to  region  can  often  be  easily 
detected. 


The  procedure  for  employing  the  multi-mode  method  begins  by  obtain¬ 
ing  spectral  amplitudes  of  the  fundamental-  and  combined  higher-mode 
Rayleigh  waves  using  the  multi pie- filter  method  (Oziewonski  al.  , 
1969).  The  fundamental  mode  is  usually  delineated  without  difficulty. 
Local  amplitude  maxima  are  taken  to  correspond  to  higher  modes  if  they 
are  within  ±0.1  km/sec  of  the  theoretical  group  velocity  for  the 
appropriate  velocity  model.  These  observed  spectral  amplitudes  are  then 


compared  to  theroetical  spectra  computed  for  various 


models  utlizing 


the  formulation  of  Levshin  and  Yanson  (1971)  to  compute  spectra  and  the 


equations  of  Anderson  et  al.  (1965)  to  include  Q  effects. 


As  shown  in  the  following  section,  the  multi-mode  method  is  pri¬ 
marily  sensitive  to  Qp  ^e  upper  crust.  Therefore  that  is  the  only 
crustal  parameter  which  is  searched  for,  and  it  is  not  difficult  to 
obtain  a  satisfactory  model  by  trial  and  error.  Qur  initial  procedure, 
in  this  study,  is  to  find  the  model  which  provides  the  best  trial-and- 
error  fits  to  all  of  the  spectra.  Our  only  adjustable  parameters  are 
therefore  Qp  the  upper  crust  and  the  seismic  moment  for  each  of  the 
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events  used  in  the  study.  The  seismic  moment  is  only  a  scale  factor 
which  sets  the  level  of  the  spectra,  whereas  the  Q  model  affects  both 
the  shape  of  the  spectra  and  the  ratio  between  higher-  and  fundamental¬ 
mode  spectral  values. 

The  possible  sources  of  error  inherent  in  the  multi-mode  method  may 
be  separated  into  two  groups:  those  which  are  associated  with  uncertain¬ 
ties  in  the  source  spectrum  and  those  which  may  occur  because  of  propa¬ 
gation  effects  between  the  source  and  receiver.  Errors  due  to  source 
effects  include  uncertainties  in  focal  depth,  strike,  dip,  and  slip  of 
the  fault  on  which  an  earthquake  occurs,  as  well  as  the  effects  of  fin¬ 
iteness  of  the  fault  and  deviations  from  the  assumed  source  time  func¬ 
tion.  For  an  ideal  explosion,  the  factors  pertaining  to  the  geometry  of 
the  fault  plane  (strike,  dip,  slip,  finiteness)  do  not,  of  course,  enter 
in.  If,  however,  an  explosion  causes  motion  on  a  fault,  then  the  possi¬ 
ble  effect  of  those  factors  on  the  spectra  cannot  be  ignored.  Errors 
due  to  propagation  effects  include  those  due  to  an  incorrect  velocity 
model  and  lateral  variations  of  elastic  and  anelastic  properties  along 
the  path  of  propagation. 

Tsai  and  Aki  (1970)  showed  that  the  effects  of  uncertainties  in  the 
angles  of  strike,  dip  and  slip  upon  the  shape  of  amplitude  spectra  are 
relatively  small,  if  those  angles  are  known  to  within  15°.  The  effect 
of  uncertainties  in  focal  depth  is  larger  and  variations  in  focal  depth 
can  produce  significant  changes  in  spectral  shape.  In  the  present 
study,  the  focal  depths  are  considered  to  be  known  quite  well,  but  if 
they  are  not,  it  is  still  possible  using  the  locations  of  spectral 
minima  in  some  cases,  to  estimate  the  depth. 
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The  effect  of  source  finiteness  should  be  considered  when  using 
earthquake  sources  and  also  when  using  explosion  sources  which  produce 
significant  fault  movement.  Ben-Menahem  (1961)  showed  that  the  ampli¬ 
tude  of  the  finiteness  factor  can  be  written  as  sinX/X  ,  where  X  = 
ui  L/2(  1/v-cos  6/e).  L  and  v  are  the  length  and  rupture  velocity  of  the 
fault,  respectively,  c  is  the  phase  velocity,  uj  is  angular  frequency, 
and  is  the  angle  measured  from  the  propagation  direction.  Assuming  a 
rupture  velocity  of  3.0  km/sec  and  a  fault  length  of  1  km,  values 
representative  of  an  intraplate  earthquake  of  magnitude  4.5  (Nuttli, 
1981),  we  find  the  finiteness  effect  to  be  insignificant  at  the  periods 
used  in  this  study.  If  the  fault  is  assumed  to  be  3  km  long,  then  the 
finiteness  factor  at  the  shortest  period  used  in  this  study  (3.5  sec)  is 
about  0.64  for  the  fundamental  mode  and  0.73  for  the  first  higher  mode. 
Considering  the  uncertainty  in  our  amplitude  data,  those  factors  are  of 
little  consequence  to  our  results.  For  explosions,  there  will  be  no 
effect  of  source  finiteness  unless  the  explosion  causes  significant 
fault  movement.  Although  that  possibility  may  occur,  fault  movement 
cannot  be  very  large  in  the  cases  we  have  considered  since  little  Love 
wave  particle  motion  was  found  on  the  horizontal  component  seismograms. 

We  assume  the  source  time  function  for  earthquakes  to  be  a  step. 
That  assumption  is  consistent  with  observations  from  numerous  studies  of 
small  shallow-focus  earthquakes.  Departures  from  a  step-function 
source,  for  the  present  purpose  would  produce  differences  in  spectral 
shape  but  would  not  affect  the  ratio  of  higher-  to  fundamental-mode 
amplitude  spectra.  The  source  spectrum  for  explosions  differs  from  that 
of  earthquakes  by  including  a  small  maximum  near  the  short-period  end  of 
the  spectrum.  The  effect  of  that  maximum  on  observed  spectra  is 


discussed  in  later  paragraphs. 

The  effects  of  changes  in  the  velocity  model  upon  spectral  ampli¬ 
tudes  has  been  studied  previously  (e.g.  Mitchell,  1975).  These  effects 
have  been  found  to  be  insignificant  compared  to  those  produced  by 
changes  in  the  Q  model  as  long  as  the  velocity  model  is  reasonably  close 
to  the  true  structure.  In  the  present  study  the  velocity  models  explain 
the  surface  wave  dispersion  data  very  well,  so  uncertainties  in  the 
velocity  model  are  not  expected  to  adversely  affect  our  results.  The 
effects  of  lateral  changes  in  elastic  and  anelastic  properties  are  more 
difficult  to  assess.  Lateral  changes  may  produce  lateral  refraction, 
multi-pathlng  and  scattering  which  can  distort  the  recorded  waveforms 
and  spectra.  The  multiple-filter  method  at  least  partially  mitigates 
adverse  effects  of  those  factors  by  separating  late-arriving  phases. 
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Q_  Models 

P 

The  method  was  first  applied  to  three  regions  of  North  America:  the 
eastern  United  States,  the  Colorado  Plateau,  and  the  Basin  and  Range 
Province,  as  well  as  to  Iran  and  the  Barents  shelf.  Q  values  in  the 

P 

upper  crust  were  found  to  be  about  250  for  the  eastern  United  States, 
150  for  the  Colorado  Plateau,  and  less  than  80  for  the  Basin  and  Range, 
using  a  trial-and-error  process.  An  automatic  inversion  procedure  pro¬ 
duced  slightly  higher  values  (Figure  2).  The  value  for  the  eastern 
United  States  is  consistent  with  that  found  independently  by  other 
methods  (Herrmann  and  Mitchell,  1975). 


Upper  crustal  Qa  values  for  the  region  of  Iran  and  Turkey  were 

P 

found  to  be  similar  to  those  of  the  Basin  and  Range.  The  data  in  that 
region  required  smaller  values  in  the  lower  crust  however.  The  model 
which  provides  the  best  fit  to  amplitude  spectra  throughout  most  of  the 
Iranian  Plateau  and  Turkey  consists  of  a  13  km  thick  upper  crust  with  a 
Q_  value  of  85  overlying  a  lower  crust  with  a  Q_  value  of  300  (Figure 

P  P 

3).  Data  for  the  region  around  Meshed,  Iran,  however,  are  better 
explained  by  a  model  in  which  the  upper  crust  is  18  km  thick  and  has  an 

average  Q„  value  of  75.  The  complexity  of  several  of  the  spectral 

P 

shapes  in  this  area  indicates  substantial  lateral  complexity  of  the  cru¬ 
stal  structure  there. 


The  multi-mode  method  was  also  used  to  3tudy  the  attenuation  pro¬ 
perties  of  the  crust  throughout  the  Barents  shelf.  The  sensitivity  of 
the  spectra  to  various  source  factors  and  to  propagation  factors  which 
might  characterize  a  shelf  region  were  thoroughly  studied.  For  plane¬ 
layered  models,  the  method  was  found  to  be  very  sensitive  to  Q.  in  the 

P 
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Figure  2.  Crustal  two-layer  models  obtained  for  the  eastern  United  States,  Colorado  Plateau  and 
Basin-and-Range.  Solid  lines  indicate  models  obtained  by  formal  inversion  of  multi-mode  spectra 
and  dashed  lines  indicate  models  obtained  by  trial-and-error . 


upper  crust  and  to  Q_  of  sedimentary  layers  if  those  layers  are  thick 

P 

and  highly  attenuating.  It  is  far  less  sensitive  to  Q_  in  the  lower 

P 

crust  and  to  Q^.  The  spectra,  for  shallow  strike-slip  earthquakes,  were 
very  Insensitive  to  changes  in  strike,  slip,  and  dip  angles  of  the  fault 
throughout  the  entire  period  range  of  interest,  but  are  very  sensitive 
to  variations  in  focal  depth  except  in  the  long-period  range  of  the  fun¬ 
damental  mode  and  the  short- period  range  of  the  higher  modes.  The 
fundamental-mode,  but  not  the  higher  mode,  spectra  for  one  path  across 
the  Barents  shelf  were,  apparently,  strongly  distorted  at  short  periods 
by  lateral  heterogeneities  along  the  travel  path.  Simple  two-layer, 
trial-and-error  inversions  yield  a  model  of  the  Barents  shelf  with  an 
upper  crustal  Q_  value  of  about  80  (Figure  4).  Three- layer  inversions 

P 

with  a  sedimentary  layer,  however,  yield  an  upper  crustal  Qn  valule  of 

P 

about  250,  if  the  sediments  1  km  thick  are  characterized  by  Q  values  of 

P 


40  (Figure  5) 


water  layer  0.3  km  thick. 
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Frequency  Dependence  of  Q 
In  the  Continental  Crust 


P 


Most  studies  which  obtain  models  of  shear  wave  Q  (Op)  as  a  function 
of  depth  assume  that  that  quantity  is  independent  of  frequency.  That 
assumption  is  computationally  convenient  and  is  based  largely  upon 
laboratory  Q  measurements  made  at  low  pressures  and  temperatures.  Dur¬ 
ing  the  period  of  this  contract,  however,  evidence  has  been  found  that 

Qp  in  the  upper  crust  varies  with  frequency,  at  least  at  frequencies 
near  1  Hz.  This  result  has  implications  for  observations  of  regional 
phases  and  yield  determinations  from  them  since  it  means  lower  attenua¬ 
tion  than  that  which  would  be  predicted  by  frequency-  ndependent  models 
obtained  from  the  analysis  of  lower- frequency  surface  waves. 

The  frequency  dependence  of  Qp  f^e  crust  of  the  eastern  United 
States  has  been  studied  using  both  frequency-domain  and  time-domain 
methods.  Only  the  results  of  the  frequency-domain  method  will  be  dis¬ 
cussed  in  this  section.  The  results  of  the  time-domain  method  are  dis¬ 
cussed  in  the  following  section  on  the  synthesis  of  short-period  Lg 
waves.  It  is  shown  there  that  approximations  made  in  the  frequency- 
domain  method  (mainly  the  use  of  only  a  few  modes)  lead  to  correct  esti¬ 
mates  of  1  Hz  Lg  attenuation  coefficient  values,  but  to  underestimates 


of  absolute  amplitudes 
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Frequency-Domain  Method 

The  Rayleigh  wave  attenuation  data  of  Herrmann  and  Mitchell  (1975) 
were  inverted  to  obtain  Q  "1  models  of  the  crust  assuming  various  types 

P 

of  frequency  dependence.  The  attenuation  coefficient  values  predicted 
for  higher  modes  for  those  various  models  were  then  compared  with 
attenuation  coefficients  observed  for  1  Hz  Lg  waves  reported  by  various 
workers  in  the  eastern  and  central  United  States. 

In  computations  of  higher  mode  attenuation,  it  was  assumed  that 
observed  attenuation  coefficient  values  at  periods  of  3  seconds  and 
larger  consist  of  only  the  first  higher  mode,  that  the  value  at  2 
seconds  is  the  average  of  the  first  and  second  higher  mode  values,  and 
that  at  a  period  of  1  second,  the  observed  value  is  the  average  of  the 
attenuation  coefficient  values  for  the  first  four  higher  modes.  Since 
differences  in  the  attenuation  coefficient  values  among  the  various 
higher  modes  are  not  great  for  the  eastern  United  States,  the  results  of 
this  study  do  not  depend  critically  upon  the  number  of  modes  employed  at 
any  period. 

The  ability  to  separate  the  effects  of  frequency-dependence  of  Qn 

P 

from  that  of  depth-dependence  is  possible  because  higher  modes  sample  a 
different  depth  interval  in  the  crust  than  does  the  fundamental  mode  for 
a  given  period.  The  higher  modes  at  1  second  period  sample  the  same 
depth  interval  as  that  sampled  by  the  fundamental  mode  at  a  period  of  6 
seconds  (about  20  km).  Likewise,  the  first  higher  mode  at  a  period  of  6 
seconds  samples  a  depth  interval  similar  to  that  sampled  by  the  funda¬ 
mental  mode  at  a  period  of  25  seconds  (about  60  km). 
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The  initial  step  in  the  procedure  employed  is  to  invert  the  funda¬ 
mental  mode  data  to  obtain  models  of  as  a  function  of  depth.  If  it 

P 

is  assumed  that  compressional  wave  Q  values,  are  twice  as  large  as 

Qp,  then  the  Rayleigh  waves  attenuation  coefficient  values,  Xr,  are 

related  to  Q  "1  and  the  elastic  properties  by 
P 

(1) 

where  4  and  p  are  compressional  and  shear  wave  velocities,  Cp  is  Ray¬ 
leigh  wave  phase  plus  velocity,  T  is  period,  and  is  a  layer  index 
(adapted  from  Anderson  et  al..  1965).  The  subscripts  ui,  4,  p,  and 
refer  to  frequency,  compressional  velocity,  shear  velocity,  or  density 
being  held  constant.  The  assumption  that  _  2Qp  has  little  affect  on 
our  results  since  surface  wave  attenuation  is  so  slightly  affected  by 

The  above  equation  is  usually  employed  assuming  that  Qp  is  indepen¬ 
dent  of  frequency.  In  the  present  study,  however,  it  will  be  assumed 

that  Qp  depend  upon  both  depth  and  frequency,  so  that  in  each  layer,  , 

(2) 

where  ^  can  be  constant  through  the  entire  frequency  range  of  interest 
or  can  vary  with  frequency.  Mote  that  if  ^  =  0,  then  Qp  is  independent 
of  frequency  in  the  commonly  assumed  way. 

The  fundamental-mode  data  are  inverted  assuming  various  values  for 
The  fits  of  the  theoretical  curves  to  the  data  of  Figure  1  for  ^  = 
0.0  and  0.5  indicate  that  the  fundamental  mode  data  can  be  explained 
equally  well  by  a  wide  range  of  modes  corresponding  to  different 
values.  This  situation  occurs  because  variations  in  frequency  depen¬ 
dence  represented  by  different  values  of  $  in  (1)  can  be  offset  by  vari¬ 
ations  in  values  of  C  at  different  depths  in  the  model. 


k 
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The  models  obtained  can  now  be  used  to  compute  attenuation  coeffi¬ 
cient  values  for  the  higher  modes  with  a  forward  calculation  using  equa¬ 
tion  (1).  The  resulting  values  can  then  be  compared  with  the  higher 
mode  data  of  Figure  1 . 


Attenuation  Coefficient  (XiOTs),  km*1 
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Results 

The  results  of  computations  for  five  values  of  between  0.0  and 
1.0,  appear  in  Figure  1.  It  is  immediately  clear  that  a  frequency- 
independent  model  ($  =  0.0)  does  not  even  come  close  to  explaining  the 
1-second  Lg  data.  Although  ^  values  between  0.3  and  0.5  provide  better 
fits,  none  of  these  curves  for  which  ^  is  constant  over  the  entire  fre¬ 
quency  range  provide  very  good  fits  to  the  entire  data  set  between  1  and 
10  seconds. 

It  is  interesting  to  note  that  all  curves  for  ^  between  0.0  and  0.5 

are  consistent  with  the  data  at  periods  of  4  seconds  and  more.  This 

consistency  supports  the  frequency-independent  internal  friction  models 
of  Mitchell  (1973)  and  Herrmann  and  Mitchell  (1975)  as  being  representa¬ 
tive  of  the  crust  in  eastern  North  America.  It  can  also  be  seen  that 
small  negative  values  of  ^  would  also  explain  the  data  at  periods 
greater  than  4  seconds. 

Further  calculations  were  performed  in  which  ^  was  allowed  to  vary 
with  period.  Although  it  is  possible  to  achieve  excellent  fits  to  the 
data  in  some  of  these  cases,  it  was  considered  that  because  of  the  large 

uncertainties  in  the  data,  that  it  is  sufficient  to  consider  only  cases 

where  ^  is  oonstant  in  this  report. 

Models  for  which  is  constant  over  the  entire  period  range  appear 
in  Figure  2.  Although  it  does  not  adequately  explain  the  higher-mode 
data,  a  frequency-independent  model  (c(  =  0.0)  is  presented,  along  with 
its  standard  deviations.  A  second  model  appears  for  the  case  where 
assumes  the  constant  value  of  0.3  at  all  periods.  Other  models  are  pos- 
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slble,  and  perhaps  even  preferable,  In  which  ^  assumes  other  constant 
values  or  varies  with  frequency.  In  the  following  section  on  the  syn¬ 
thesis  of  short-period  Lg,  it  appears  that  y  =  0.5  explains  synthetic 
data  reasonably  well. 

The  computations  of  the  present  study  indicate  that  in  the 

P 

crust  in  eastern  North  America  is  frequency  dependent.  Several  possible 
models  can,  however,  explain  fundamental-  and  higher-mode  Rayleigh  wave 
attenuation  coefficient  data  at  periods  of  1  second  and  greater.  The 
models  differ  in  the  nature  of  their  dependence  upon  frequency.  One 
reason  for  the  non-uniqueness  of  the  results  is  the  occurrence  of  a  gap 
in  the  higher-mode  data  at  periods  between  1  and  4  seconds,  if  reliable 
data  can  be  obtained  for  that  range,  it  would  be  possible  to  exclude 
some  of  the  models  we  have  considered.  Later  work  (Mitchell,  1981)  does 
indeed  suggest  that  we  can  eliminate  models  which  produce  peaks  in 
attenuation  at  periods  between  1  and  4  seconds. 
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Synthesis  of  Short-period  Lg 
in  the  eastern  United  States 
by 

Brian  J.  Mitchell 

Recent  work  (e.g.  Bache  si.  al. ,  1981;  Wang,  1981)  has  allowed  us, 
for  the  first  time,  to  compute  short-period  seismograms  utilizing  modal 
summation.  This  allows  us  to  compare  computed  Lg  waveforms  with  those 
actually  observed  on  short- period  records. 

Early  computations  were  performed  to  see  how  many  modes  are 
required  to  obtain  useful  attenuation  coefficient  values.  It  was  found 
for  sources  more  than  a  few  km  deep  that  as  few  as  5  modes  gave  useful 
values  for  an  eastern  United  States  model.  For  very  shallow  source 
depths  (less  than  1  km)  and  a  thick  low- velocity  sediment  layer,  more 
modes  are  needed,  but  10  modes  seem  sufficient  In  all  cases.  Absolute 
amplitude  values,  however,  require  more  modes  to  be  realistic.  It 
appears  that  at  least  15  modes  are  needed  and,  in  some  cases  at  least, 
as  many  as  30  or  more  are  desirable  to  obtain  realistic  results. 

Later  work  showed  that  it  was  necessary  to  include  a  surface  layer 
with  relatively  low  velocities  if  we  are  to  compute  realistic  short- 
period  seismograms  of  the  Lg  phase.  Figures  1  and  2  compare  a  model 
having  low-velocity  surface  sediments  with  a  model  lacking  such  a  layer 
(see  Table  1).  The  latter  model  produces  seismograms  having  spikelike 
arrivals  which  are  not  observed  on  real  seismograms. 


Figure  2  •  Vertical  component  synthetic  seismograms  at  4  distances 
for  the  velocity  model  of  McEvilly  (1964).  Qa  ■  2000  and 
Qj3  ■  1000  at  all  depths. 


68 


Tabie  1 


Eastern  and  Central  United  States  Velocity  Models 


Simplified  EUS  Model 


h 

4 

B 

1 

P 

C  .25 

2.42 

1 .30 

2.30 

0.75 

it  .91 

2.71 

2.6? 

6  .00 

8.19 

2.4? 

2.75 

18.  CO 

6  .38 

3.67 

2 .90 

18. 0C 

7.17 

3.98 

3.00 

— 

8. 1C 

'1.58 

3.30 

McEvil ly 

(  196*0 

h 

4 

P 

P 

1 1  .00 

6 .10 

3.50 

2.70 

9.00 

6  .40 

3.68 

2.90 

18.CC 

6  .70 

3.67 

2.90 

24. CO 

8.15 

4.67 

3.30 

^C.OO 

8.20 

4.47 

3.30 

16C.00 

8.20 

4.45 

3.30 

— 

8.70 

4.80 

3.60 
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Comparison  with  Observed  Seismograms 

Seismograms  recorded  at  several  stations  for  a  small  central  United 
States  earthquake  of  14  August  1965  were  digitized  for  comparison  with 
synthetic  records.  Synthetics  were  computed  for  three  models:  one  in 

which  0^  ia  independent  of  frequency,  one  in  which  Qp  varies  with  fre¬ 
quency  as  uP'5,  and  one  in  which  varies  vith  frequency  as  ui0,5  but 
which  includes  a  layer  of  low-Q  sediments.  Examples  are  shown  in  Fig¬ 
ures  3  and  4  for  Rayleigh  waves  and  Love  waves,  respectively.  It  can  be 
seen  that  most,  if  not  all,  of  the  observed  coda  can  be  explained  by 
synthetic  seismograms  in  plane-layered  models.  It  is  possible  to 
explain  the  gross  features  of  the  coda,  such  as  duration,  absolute 
amplitude,  and  perhaps  envelope  shape  in  some  cases,  but  not  the  details 
of  the  wave  form. 

At  short  distances,  all  frequency-dependent  models  tested  produced 
a  more  prominent  and  more  regular  fundamental-mode  Rayleigh  wave  train 
than  that  which  was  observed.  At  larger  distances  the  time  interval  on 
observed  records  corresponding  to  the  fundamental  mode  is  comprised  of 
more  irregular  wave  forms  than  the  synthetics.  This  time  interval  may 
contain  the  remnant  of  the  fundamental  mode  after  attenuation  by  intrin¬ 
sic  Q  and  scattering  as  well  as  converted  modes  -  fundamental  to  higher 
modes  -  and  vice  versa. 

Figures  5  and  6  present  comparison  of  observed  and  theoretical 
amplitudes  from  Rayleigh  and  Love  waves,  respectively,  for  two  Q  models. 
A  frequenoy-dependent  model  in  which  Q  varies  with  frequency  as  uP»5 
provides  a  better  match  in  both  cases. 
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Figure  5.  Observed  Rayleigh  amplitudes  at  1  Hz  versus  theoretical 
amplitudes  predicted  by  frequency-independent  (5  ■  0-0)  and 
frequency-dependent  (5  -  0.5)  Qg  models  of  Mitchell  (1980). 


Figure  6.  Observed  Love  amplitui 
predicted  by  frequency-indej 
(1980). 
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